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ABSTRACT 

Recent studies of heavy r-process elements in low [Fe/H] halo stars have suggested that an initial 
population of metal-free very massive stars (VMSs) may be required to provide early Fe enrichment 
without coproducing heavy r nuclei. We find similar abundance trends in a-elements (which should be 
copiously produced by VMSs), but not in other elements such as carbon (which should not), in agreement 
with this hypothesis. We then combine the corresponding level of prompt initial enrichment with models 
of VMS nucleosynthetic yields and spectra to estimate the corresponding ionizing fluxes. The result 
suggests that there may have been enough VMS activity to reionize the universe. The unusually hard 
spectrum of VMSs would imply a different reionization history from canonical models. Hen could have 
been reionized at high redshift, only to recombine as a subsequent generation of stars formed with a 
"normal" initial mass function. 

Subject headings: cosmology: theory - early universe - galaxies: formation -stars: abundances 



1. INTRODUCTION 

The nature of the very first generation of metal-free 
('Population IIP) stars has been debated for decades 
(Schwarzschild & Spitzer 1953; Ezer & Cameron 1971). 
It is often thought that the lack of efficient cooling mech- 
anisms would lead to a top-heavy initial mass function 
(IMF) (Larson 1998), and particularly to the production 
of 'very massive stars' (VMSs) with M* > 100 Mq (Carr, 
Bond & Arnett 1984). Recent numerical simulations of 
cooling primordial molecular gas have hinted that this 
is indeed the case: they find quasistatically-contracting 
gas clumps with characteristic densities ^ 10* cm"-^ and 
temperatures of a few hundred kelvin - a regime set by 
the properties of the chief coolant, molecular hydrogen 
(Bromm, Coppi, & Larson 2001a; Abel, Bryan & Nor- 
man 2000; Nakamura & Umemura 2001). The Jeans un- 
stable clumps are of order a few hundred solar masses; 
because the contraction is subsonic, further fragmentation 
is unlikely, suggesting that the very first stars were indeed 
very massive. Among other properties, VMSs would pro- 
duce many more ionizing photons per unit stellar mass 
than ordinary stars (Bromm, Kudritzki, & Loeb 2001c), 
and have much harder spectra; they could thus ionize Hen 
(Tumlinson & ShuU 2000) . Their existence would have im- 
portant implications for the metal enrichment and reion- 
ization history of the universe. 

The duration and extent of early metal-free star forma- 
tion are not known, and depend on the complex and in- 
homogeneous process of early metal enrichment (Gncdin 
& Ostriker 1997; Madau, Ferrara, & Rees 2001). At 
present, the best evidence for VMSs with M* < 250 Mq 
may be from their distinctive nucleosynthetic products: 
VMSs process more of their mass into alpha-elements than 
normal stars, and they finally explode by a mechanism 
(pair instability) that is very different from core-collapse 
supernovae (Fryer et al. 2001), and which is unlikely to 
produce r-process elements (Heger & Woosley 2001, here- 



after HW) . Metal-poor halo stars provide important fossil 
evidence for early nucleosynthesis because they extend to 
metallicities even lower than those inferred for the damped 
Lya systems (e.g., Ryan 2000). 

Recent observations of r-process abundances as a func- 
tion of [Fe/H] in halo stars suggest that an initial pop- 
ulation of rapidly evolving VMSs provided an intial Fe 
inventory without co-producing heavy r-nuclei (Wasser- 
burg & Qian 2000, hereafter WQ; Qian & Wasserburg 
2001b). In this Letter, we estimate the total number of 
VMSs necessary to produce the required pre-enrichment, 
and consider the consequences for reionization, in the 
spirit of previous papers (e.g. Haiman & Loeb 1997, 
Miralda-Escude & Rees 1998) which have considered the 
connection between the metallicity of the Lya forest 
at z ^ 3 and high-redshift ionizing photon production 
for a 'normal' IMF. In all numerical estimates, we as- 
sume a ACDM cosmology with (ilM, f^A, ^b, h, cr^f^-i , n) = 
(0.35, 0.65, 0.04, 0.65, 0.87, 0.96). 

2. ENRICHMENT BY VERY MASSIVE STARS 

In Figure |^ we show VMS yields as functions of mass, 
as computed by HW for stars in the mass range subject to 
pair instability explosions Af* 140-260 Mq. Such stars 
are completely disrupted and leave no remnants. Stars 
with 40 M© < M* < 100 Mq and > 250 Mq are pre- 
dicted instead to collapse to black holes without metal 
ejection (HW); thus, nucleosynthetic tracers only provide 
a lower limit on early star formation. The top panel shows 
the yields yi = Mi/M^, (where Mi is the total mass of 
element i produced, and is the initial stellar mass), 
assuming a He core mass of Mcoro ~ 0.5 as in HW. 
In zero metallicity stars the efficiency of nuclear-powered 
radial pulsations (Baraffe, Heger, & Woosley 2001) and 
radiatively driven stellar winds (Kudritzki 2000) is greatly 
reduced, so mass loss is likely to be small. From the figure 
it is clear that iron provides an uncertain measure of VMS 



^Theoretical Astrophysics, Mail Code 130-33, Caltech, Pasadena, CA 91125 

^Department of Astronomy and Astrophysics, University of California Santa Cruz, Santa Cruz, CA 95064 
^The Lunatic Asylum, Division of Geological and Planetary Sciences, Caltech, Pasadena, CA 91125 

1 



2 



Did Very Massive Stars Pre-enrich and Reionize the Universe? 



formation. The amount of silicon burned into iron in the 
final VMS explosion increases with the kinetic energy of 
the initial collapse, and thus the iron yield varies strongly 
with the core mass. Yields of other nuclides depend less 
on M*. The lower panel of Fig. |^ depicts the ratio of 
the mass of species i produced in the VMS to the total 
mass of metals, divided by the analogous ratio for solar 
abundances. It illustrates the non-solar yields of VMSs. 
While the mass yield of oxygen is quite high, the a ele- 
ments Si, S, and Ca have enrichment factors relative to 
solar which are far above O. These trends are specific to 
very low/zero-metallicity VMSs. The Si yield is about 
- 10% in a - 150 M© VMS which is metal-free (HW), or 
has a metallicity Z = 4 x IO^^Zq (Portinari et al. 1998). 
However, the Si yield drops by between 2 and 4 orders of 
magnitude for VMSs with metallicity Z > AxlO'^ Zq (Ta- 
ble 12 of Portinari et al. 1998); likewise, the maximum Si 
yield of Type II supernovae (SNe) from progenitor stars of 
any mass and metallicity is about an order of magnitude 
less, 1%. Thus, at low metallicities Si, S, Ca abun- 
dances should be excellent tracers of Population III VMSs. 
No elements above the iron group are produced in the HW 
models for the pair instability explosion. If these models 
are correct, VMSs should produce copious a-elements with 
very little associated r-element production. 
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Fig. 1. — VMS yields from Heger & Woosley (2001) used in 
this paper as a function of the He core mass. Top panel: the yield 
yi = Mi/Mt, , where Mi is the mass of species i ejected, and M* is the 
initial stellar mass [roughly twice the final He core mass M(core)]. 
Bottom panel: See text. 

We now turn to observations. WQ analysed data from 
observations of heavy r~process abundances in metal-poor 
halo stars as a function of [Fe/H] (see Fig. ||, top left). 
They found the following trend in the data. At the lowest 
metallicities, [Fe/H]< —3, the Ba abundance is roughly 
constant. Then, over a small range in Fe abundance, 
—3.1 < [Fe/H] < 2.5, there is a wide spread in the Ba 
abundance of order 2 dex. They interpreted this sudden in- 
crease in the Ba abundance dispersion as follows. Previous 



studies of radioactive ^^^I and ^*^Hf in meteorites appear 
to require two distinct classes of Type II SNe to account for 
r-process elemental abundances: high-frequency events 
(H) with characteristic timescales ~ lO^yrs that produce 
heavy r-process elements with mass numbers A > 130 
(e.g. Ba, Eu) and low-frequency events (L) with character- 
istic timescales ~ lO^yrs that produce light r-process ele- 
ments with mass numbers A < 130 (e.g. iodine) (Wasser- 
burg, Busso, & Gallino 1996). Because of their shorter life- 
time, the H events are the first SNe to explode. As they 
are not associated with significant coproduction of iron, 
they will produce a wide scatter in Ba abundance over 
a small range in [Fe/H]. However, the region [Fe/H]< —3, 
where the Ba abundance is very low and roughly constant, 
was considered by WQ and Qian & Wasserburg (2001b) to 
be an initial or prompt inventory dominated by the con- 
tributions of very low-metallicity massive first-generation 
stars, which appear before H events and produce very 
small amounts of r-process elements relative to their iron 
production rate. In this case, the iron enrichment up to 
[Fe/H]'^ —3 is almost exclusively due to VMSs. The finite 
Ba abundance at [Fe/H]< —3 can be associated with very 
minor coproduction of Ba in VMSs, or trace contamination 
by H events (see WQ and Qian & Wasserburg 2001b). 
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Fig. 2. — Observed Ba abundances vs. other metal abundances 
in metal poor halo stars, from the data of McWilliam et al. (1995). 
The original trend found by Wasserburg & Qian (2000) was in Ba vs. 
[Fe/H]; the sharp increase in scatter of Ba/H at fixed [Fe/H]~ —3 
also occurs in elements overproduced by VMSs, such as Si and Ca, 
suggesting that the sources of Ba and Fe, Si, Ca elements are dis- 
tinct. Elements not strongly overproduced by VMSs, such as C, do 
not show such a trend. 

From our previous discussion we would expect the trend 
seen in Fe to be most strongly exhibited in the a elements. 
In Fig.||, we plot halo-star Ba/H values against a element 
abundances, using the data of McWilliam et al. (1995). 
Indeed, Si and Ca show the same trend of a low Ba abun- 
dance at very low metallicity followed by a wide scatter 
over a small range in metallicity. On the other hand, we 
would expect these trends to be absent in elements less 
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abundantly produced by VMSs, which do not experience 
significant initial enrichment before normal star formation 
takes over. This is indeed seen: carbon shows merely a 
linear correlation between C and Ba production with a 
great deal of scatter (lower right panel). This provides 
additional supporting evidence for early VMS formation. 
Most importantly, while VMS Fe yields vary strongly as a 
function of mass. Si and Ca yields do not change widely 
above Mcnrc > TOM© (M* > 150 M©, from from eqn. 1 of 
HW) . In §11 we will use these yields to quantify the amount 
of early VMS formation. 

Although oxygen is the element produced most abun- 
dantly in VMSs, unfortunately we cannot use it to mea- 
sure VMS formation. Oxygen abundance estimates in low 
metallicity stars are uncertain, with different trends of 
[0/Fe] vs. [Fe/H] claimed by different observers (e.g., Is- 
raelian et al. 2001; Carretta et al. 2000). Moreover, to 
the best of our knowledge, there is no substantial body of 
simultaneous O and r-process abundances in halo stars. 
Note that the predicted O yields of VMSs are similar to 
SNII yields, and so the nucleosynthetic signature of VMSs 
may be less distinct in O than in the a elements (see, 
however, Qian & Wasserburg 2001a). 

3. CAN VERY MASSIVE STARS REIONIZE THE UNIVERSE? 

VMSs are ~l-2 orders of magnitude more efficient than 
normal stars with a standard IMF at producing ioniz- 
ing photons, per unit stellar mass (Bromm et al 2001c). 
Thus, even if the cooling and star formation efficiencies in 
gas with primordial composition are low, a much smaller 
fraction of the gas needs to be turned into stars in or- 
der to reionize the universe. We now attempt to deter- 
mine if VMSs could indeed have reionized the universe. 
The onset of the large increase in loge(Ba) occurs when 
the abundances of Si and Ca are [Si/Hjjump w —2.3 and 
[Ca/Hjjump ~ —2.4. Let us assume that the IGM was uni- 
formly polluted by VMSs up to this critical metallicity. (Si 
abundances in the Ly-a forest yield [Si/H] ~ —2.3, with 
a scatter of a factor ~ 3 [Cowie & Songalia, 1998; Elli- 
son et al., 2000], in good agreement with our estimate of 
[Si/HjjumpO The number density rii of species i relative 
to hydrogen {ni/nn) at the critical metallicity is deter- 
mined by the fraction /vms of baryonic matter that was 

processed in VMSs. Using [i/i^ljump = log (^T^^^yf^jf^) , 

we have 

/vms = (rii /nn) jump /{ni/nH)YMS = ^i,o x lOl'/^lj"'"''/?/^ 

Here, Xi^Q is the mass fraction of species i in the sun 
and Ui is the mass yield of species i from a VMS. From 
Fig. 1, we have ygi — 0.05 — 0.1 for all the stars with 
15OM0 < M, < 25OM0. We thus estimate that a frac- 
tion /vms ~ 3 — 7 X 10^^ of all baryons must have been 
processed through VMSs. Ca yields a similar range, but 
extending up to 1.2 x 10^**. The result is fairly insensitive 
to the assumed IMF since Ca and in particular Si yields 
are fairly constant for Mcorc > 70Mq (M* > I5OM0). 
An IMF dominated by VMSs with < I5OM0 can be 
excluded since these stars produce very little Fe, and the 
observed jump in Ba abundance at [Fe/II]~ —3 would not 
be seen. 

In the above estimate we assume that metals produced 
by VMSs are uniformly mixed throughout the IGM, and 



this is a fairly uncertain assumption. It is possible that 
early star formation was highly biased and VMSs only 
polluted some small fraction of the IGM from which the 
halo stars formed. Our estimate of /vms would then be 
correspondingly reduced. Nonetheless, metals have been 
detected even in underdense gas in the IGM at 2 < z < 3, 
far from star forming regions (Schaye et al 2000), plausibly 
the result of early uniform enrichment by PopIII stars. We 
can use extended Prcss-Schechter theory to compute the 
extent to which VMS debris should be over-represented in 
the Galactic halo; for typical assumptions, the enhance- 
ment factor is only ~ 1.6 (Madau & Rees 2001). 

How many ionizing photons does our estimated /vms 
correspond to? Bromm et al (2001c) have computed mod- 
els of VMS spectra and find that for 300 M0 < M* < 
1000 M0 VMSs, the stellar effective temperature and lu- 
minosity per solar mass are approximately constant. From 
their models, a fraction /*,bkl ^ (10, 1,4) x 10^^ of the 
baryons in the universe must have been processed through 
VMSs to obtain 10 (HI,HeI,HeII) ionizing photons per H 
and He nucleus in the universe respectively (see their Ta- 
ble 1 and eq. 9). From this, we estimate that VMSs must 
have radiated /photon ^ 10/vms//*,bkl '^(3-7, 30-70, 8- 
15) photons above the (HI,HeI,HeII) ionization threshold, 
per H or He nucleus. This is probably an overestimate, 
because the stars which eject metals have lower masses, 
150 — 250 M0, than the stars considered by Bromm et al. 
(2001b). For 100 M© stars the number of ionizing pho- 
tons per solar mass is reduced by a factor ^ 1.5 — 2 for 
HI and Hel and ^ 4 for Hell (Bromm et al 2001c; Tum- 
linson & ShuU 2000). The uncertainties in our estimate 
are almost certain to be dominated by our assumption of 
uniform metal mixing, and by our neglect of more massive 
stars > 260 M0 which do not eject metals but collapse 
to black holes. 

Could these ionizing photons escape from their host 
halos? Suppose that VMSs originated in the first pre- 
galactic objects able to cool within a Hubble time, i.e. 

in halos of mass M > Merit = 4 x 10^ (^)~^ Mq, (Fuller 
& Couchman 2000; Haiman, Thoul & Loeb 1996). Gas 
within such small halos is quickly ionized by the fiux 
from a single VMS, and the escape fraction of ionizing 
photons is therefore likely to be of order unity. We ob- 
tain the density profile of the gas n{r) by assuming it 
is initially isothermal at the virial temperature and de- 
manding hydrostatic equilibrium within an NFW dark 
matter profile (Makino et al 1998). The maximum re- 
combination rate is TVrcc ~ ctB j n^{r)4TTr^dr = 3 x 
10*^ ( skms-O ^ (W)^^^s-i. This is much smaller than 
the VMS HI ionizing photon production rate, A^ion ~ 
5 X 10*^° ( 300 Mq ) photons s~^, and so the halo gas re- 
mains wholly ionized. Moreover, gas within these small 
halos is easily photoevaporated: the baryonic gravitational 

binding energy, E, 7 x lO^^ (^^^f^j'^' [(1 + z)/20] 

erg, is considerably less than the total ionizing radiation 
produced by a VMS, Eion ~ 105''(M*/200 M0)erg. The 
entire baryonic component of the halo will be photoe- 
vaporated on the sound-crossing timescale of the halo, 

icvap |rvir/c. ~ 3 X 106(53^) (i±i)-^/Vrs, com- 
parable to the lifetime of the VMS, ~ 3 x lO^r. 
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The amount of reionization is determined by competi- 
tion between ionization and recombination. The increased 
mean density of the IGM at high redshift n cx {1 + z)^ 
significantly increases recombination rates. On the other 
hand, the gas clumping factor due to dense non-linear ob- 
jects C = (riff) / {uh)^ declines very sharply at high red- 
shift, and is C* « 2 at z = 20, compared to C « 30 at 
z = 10 (Haiman, Abel & Madau 2001). Hence, the re- 
combination time ^rcc = i/o^BiT-eC is in fact longer by 
a factor - (30/2)/(20/10)3 ~ 2 at z = 20. If (many 
generations of) VMSs form over a Hubble time in at 
z ^ 20, then the number of recombinations is n^ec ^ 
tn/U-cc « (9,50)[(1 + z)/20]i-5(C/2) per Hii and Hem 
particle, respectively. Since the halo star abundances sug- 
gest that n.ion ^ (3 — 7, 30 — 70, 8 — 15) photons above 
the (HI,HeI,Hen) ionization threshold per H or He nu- 
cleus were produced, the filling factor of Hii and Hem 
regions due to VMSs may be ~ riion/nrec ~ 0.35 — 0.75 
and ~ 0.15 — 0.30 respectively. We may have underesti- 
mated the ionizing photon production rate and thus the 
filling factors by factors of a few, since we do not count 
VMSs which collapse to black holes and do not eject their 
metals. It is thus plausible that the IGM was completely 
ionized by VMSs. Once the IGM is reionized, a much lower 
comoving emissivity is required to keep it ionized at lower 
redshifts than to ionize it ab initio: the high IGM tem- 
perature strongly suppresses gas accretion onto low mass 
halos Vc < 30 km s^^, so the gas clumping factor remains 
small at lower redshift. This could reduce the comoving 
emissivity required for an ionized IGM by a factor C ^ 30 
at z = 10. 

Our estimate that /vms ~ 3 — 7 x 10^^ is consistent 
with other estimates. Due to feedback processes, it is un- 
likely that more than one VMS forms per halo (Madau 
& Rees 2001). If VMSs form in the first pregalactic 
objects able to cool within a Hubble time, with masses 
M > Merit = 4 X lO'^ (^)~^Mq, at a time when such ha- 
los collapse from 2.5 cr density peaks in the universe and 
contain ~ 1.3% of the total mass of the universe, then 
/vms ~ 1-3 X 10-2 X A/vMs/(/BMerit) = 5 X 10-5 (where 
/s = Qb/^m)- Normal star formation probably starts 



when the gas metallicity is sufficiently high for efficient 
metal-line cooling, reducing the Jeans mass and enhanc- 
ing fragmentation. Bromm et al. (2001b) find from nu- 
merical simulations that this is likely to take place once 
Z reaches 1 to 10 x IO'^Zq ~ 2 to 20 x 10"*^. Similarly, 
Hellsten & Lin (1997) find that the critical metallicity for 
isobarically cooling gas at nT ^ 10^ K cm^'^ to cool to 
T ~ 50K is [Fe/H]'^ —3, the critical metallicity identified 
by WQ. Since half the main-sequence mass of a VMS be- 
comes metals in the HW models, our estimate of fvMS 
corresponds roughly to Z ^ 0-5fvMS ~ 1-5 to 3.5 x 10"^. 
However, caution is necessary because the cooling studies 
used solar metal abundances scaled with metallicity, which 
are very different from the HW yields. The difference may 
be unimportant if oxygen is the dominant coolant (since 
O yields from HW scale with metallicity in roughly solar 
ratios, see Fig. 1), but a direct comparison requires a rig- 
orous calculation of the cooling properties of gases with 
non-solar ratios of metal abundances. 

Finally, due to the high effective temperature of VMSs, 
Hell may reionize at high redshift in VMS reionization sce- 
narios, only to recombine once ordinary stars with softer 
spectra become predominant (subsequent Hell reioniza- 
tion only took place at z ~ 3, as observed by Kriss et 
al 2001). In standard stellar reionization scenarios, Hell 
reionization takes place only once, when quasars turn on. 
The higher temperatures achievable with early Hell reion- 
ization would have interesting consequences for the ther- 
mal history of the IGM (Miralda-Escude & Rees 1994; 
Abel & Haehnelt 1999), and the Hell recombination lines 
produced in the host halo may eventually be observable 
with NGST (Oh, Haiman & Rees 2001; Tumlinson, Giroux 
& ShuU 2001). 
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